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Abstract 

The multichannel generalization of the theory of spectral, scattering and 
decay control is presented. New universal algorithms of construction of com- 
plex quantum systems with given properties are suggested. Particularly, trans- 
formations of interaction matrices leading to the concentration of waves in a 
chosen partial channel and spatial localization are shown. The limiting in- 
structive cases illustrating different phenomena which ooccur with the combi- 
nation of 'incompatible' properties are considered. For example, the scattering 
solutions with different resonance widths at the same energy for the same in- 
teraction are revealed. Analogously, a 'paradoxical' coexistence of both strong 
reflection and absolute transparency is explained. The case of the violation 
of 'natural' asymptotic behavior of partial wave function is demonstrated : 
it has a greater damping decrement for the channel with a lower threshold. 
Peculiarities of the multichannel periodic structures, bound states embedded 
into continuum, resonance tunneling and degeneracy of states are described. 



1 INTRODUCTION 

Many years ago this journal published Feshbach's paper "Unified theory of nu- 
clear reactions" [1] which attracted attention of a wide physical community. This 
approach is a mighty tool for describing quantum systems (e.g. few-body) with 
internal degrees of freedom. It can simplify the Schrodinger equation in partial 
derivatives reducing it to a system of coupled ordinary differential equations for the 
vector-valued wave functions, see (0). 

The theory has gradually been developing in the direction of its greater unifi- 
cation. Generalization to reactions with the rearrangement of particles has been 
suggested by us in this journal [2-4]. A new contribution to the theory consists 
in employing the inverse problem (IP) and quantum supersymmetry (SUSYQ) ap- 
proach. We get a deeper insight into a relationship between interactions (V) and 
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observables S = {En, Can}'- spectral weights Can or Man, energy levels and their 
analogs, e.g., resonance positions. 

In the IP and SUSYQ, the input data are S. This is a remarkable advantage. 
It provides us with a complete set of elementary quantum transformations corre- 
sponding to variations of individual spectral parameters. They are given by explicit 
expressions: the exactly solvable models, see, e.g., [5-7] and references therein for 
the one-channel case. For comparison, in the direct problem, variation of only one 
chosen S-parameter is an extremely difficult problem. 

There is in principle a possibility of converting any system into another one of a 
different nature. Specifically, a scattering state can be 'raked up' into a bound state 
embedded into continuum (S-transformation at a single energy point without any 
other spectral variation), see below. 

A lot of qualitative physical information is revealed via a computer visualization 
of these transformations. So, in the one-dimensional and one-channel case we have 
developed a qualitative theory of quantum design [8-12]. It becomes clear what 
"bricks and blocks" are needed to build systems with given physical properties. A 
comparison with a children's toy constructor-set suggests itself. In this way, the 
understanding of wave mechanics is raised to the intuitive level. We acquire the 
possibihty of a quahtative foresight of the results of IP-SUSYQ transformations 
even without formulae and computer calculations. Namely, we look for universal 
rules of the structural and dynamic response of quantum complexes when varying 
S. Particularly, an arbitrary energy level En can be created, removed or shifted. 
Any object is uniquely determined by a complete set of S data. So, it allows the 
choice of interactions providing the prescribed properties of the system. 

The same status for the multichannel theory is still to be achieved. The one- 
channel IP scalar transformations represent a zero-measure subset of their multi- 
channel matrix analogs. Thus, one can expect substantial advancement and greater 
diversity of new effects in comparison with the one-channel case. 

In this paper, the latest results on the multichannel quantum design are dis- 
cussed. First, we give short preliminary information about the systems of coupled 
Schrodinger and IP equations. 

2 Multichannel equations of direct and inverse 
problems 

The description of a wide class of quantum objects with the multidimensional wave 
functions ^' can be reduced to the system of close coupled Schrodinger (multichan- 
nel) equations if we choose one variable x and expand ^ over the set of known basis 
functions ^a{0 other variables ^ 



The coefficients of these expansions, called channel wave functions, are determined 
by the system of equations 




(1) 



a 



Ki^) + E^a/?(^)*/3(^) = Ea^aix), Ea ^ E - 6, 
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v^p{x) = J ^:{ov{x,o^pm^ (3) 

where V{x,C,) is the potential for the corresponding mult i- dimensional problem and 
Va/3{x) is the interaction matrix with non-diagonal elements coupling the equations. 
Greek indices mean the channel numbers, i.e., partial equation numbers. Sym- 
bol E stands for a total energy value and are the threshold energies. The potential 
matrices Va/3{x) depend on x locally. 

o o 

The IP gives the recipe of transformation from the initial (x), Vap (x) to the 
new Va/3{x) with the given spectral data - see, e.g., for Eqs. (|)-(0) 

oo 

=ia (X) + J J2K^i3{x,x') {x')dx'] (4) 



d 

Vapix) =Vaf3 (x) - 2—Kaf3{x, x). (5) 

Here, the transformation matrix kernel K is determined by the IP equation 

oo 

Kap{x, x') + Qapix, x')+ / ^ Ka^{x, y)Q.y(i{y, x')dy = 0. (6) 

i ^ 

The input spectral information is contained in its matrix valued kernel 

Qal3{x,x') = 

oo 

J hdhY^Fap, {x,k){k-'/'[kk)-SCk)]k-'/'}p:^ Fpa' {x',k) 



-oo 



a'/3' 



n a'l3' 



o o 



En) Ma'nMp'nFpa' {x' , En)- (7) 

In this expression the partial spectral parameter Man stands for the factor in the 
asymptotic behavior of the partial wave in the ath channel of the normalized bound 
state at the energy En- \l/a(a;, — >• Mq,„ exp(— -y/e^ — En x), as x —>■ oo. We 
use Fap{x,En) to designate the matrix Jost solutions of (0) defined at the energies 
of bound states: Fa/3{x,En) Sap exp(— A/e^ — En x), x — > oo. Here, the first 
index means the channel number (i.e. partial equation number in (§)), the second 
one designates the type of a boundary condition. In (|^ the scattering matrix is 
S(k) and Fai3{x, k) is the matrix Jost solution for the continuous spectrum. The 
sums in Q correspond to energy levels of bound states and the integral refers to the 
continuous spectrum. Symbol /c is a diagonal momentum matrix kaSa/B, ka = y/Ea- 
If initial and final systems differ in only a finite number of discrete spectral 
parameters, the kernel Q becomes degenerate and equation reduces to exactly 
solvable algebraic equations (see, e.g. fl^). Indeed, in this case the integral term 

o o o 

in M) vanishes because S = S and the terms in the sums with En = En, Man 
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= Man are cancelled. So, there remain a finite number of terms with the factorized 
dependence on x, x'. Thus, IP gives rise to an infinite number of exactly solvable 
models which form a complete set, i.e. there is a possibility (at least theoretically) 
to fit any quantum system by them. 

The above formulae correspond to the case when the spectral parameters {S, 
M) characterize the asymptotic {x oo) behavior of the wave functions. This 
is Marchenko's approach (M). It is often convenient also to control the spectral 
parameters which determine the wave behavior at the origin according to Gelfand- 
Levitan's approach (GL). Then, instead of Man the derivatives of normalized bound 
states at X = are used ^ 

Can = -J^"^ a{x, En)\x=0- 

The matrix- valued regular solutions ^a/six, E^) enter in GL formalism instead of 
the Jost ones ($„^(a;, E„) |^=o = 0, -^^apix, En)\x=o = Sa/s), and 

<iJaix, En) = ^c^pi^^ K)Cpn- (8) 

/3 

Furthermore, the minus sign in is replaced by the opposite one and the integration 
is performed over the interval [0,x]. In what follows Cj^n "will be referred to as a 
spectral weight vector (SWV) related to the nth eigenstate. 

For the reader interested in the IP approach it would be useful to conceive the IP 
as the Gram-Schmidt orthogonalization of the set of eigenvectors generalized to the 
case of infinite dimensions and continuous basis. The bound state wave functions are 
considered as vectors in a special Hilbert space in which the coordinate x numbers 
the vectors while the energy indices n indicate their components. The inner product 
in that space is determined by a measure dependent on the SWV vectors C or M 
etc, [Q. In fact, Parseval's completeness relation 

En)^>f,{x\ En) = 5a^5{x - x') (9) 

n 

can be rewritten using as follows : 

5^^$„^(x,E„)C^„Cy„$^y(x',E„) = 5ap5{x-x'). (10) 

The last equality represents the orthogonality relation with the left part as an 'inner 
product' for vectors $a^(a;,-E„). It is determined via summing over the multi-index 
(a, n) with a weight C^nCyn which gives the measure of a new Hilbert space. The 

o 

IP transformation Can^ Can means the change of the measure. Then, the Hilbert 
space is changed so that the initial orthogonality relation is violated. For this 
relation to be restored, the Gram-Schmidt orthogonalization is performed. So the 
IP equation (^ itself is just the equation for the coefficients of such a procedure. 

It turns out that the exactly solvable models of the IP can be alternatively de- 
rived within the framework of the supersymmetric quantum mechanics (SUSYQ) 
formalism by E.Witten. The SUSYQ has been generalized to the multichannel 
case [15-19]. In many respects the SUSYQ scheme is the same in both one- and mul- 
tichannel cases. The initial Hamiltonian, the second order differential operator, is 
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factorized into the simpler first order matrix operators A"^ = ±d+W{x) {d is a sym- 
bol of the derivative, ImW{x) = 0, = W{x)) [^. The transformation 
consists in permutation of : 

Ho = A+A- +S ^ Hi = A-A+ + S, (11) 

where £^ is a constant factorization energy. 

The matrix W{x) has the form as in the one-channel case ||2ll 

W{x) = %{x,S)<ilo{x,S)-' 

where \l/o(^; ^) is an arbitrary matrix solution for the initial system at the factoriza- 
tion energy. It should be noted here that the choice of £ and the selection of "ifo^x, £) 
are crucial points predetermining the resulting properties of the transformed Hamil- 
tonian and corresponding solutions. It is remarkable that there is a simple relation 
between the solutions '^q{x,E) and i?) of the Schrodinger equation with the 

initial and new Hamiltonians Hq, Hi ||22| 

E) = A-^o{x, E) = {-d + W{x))^o{x, E). (12) 

Let us perform the SUSYQ transformation twice at the same energy £ which can 
be arbitrary. The resultant formulae can be equivalent to the IP ones for variation of 
spectral weight or for creation (removal) of any bound state while all other spectral 
parameters from the complete set {i?„, C„} or M„} remain unchanged [|1^ 



see 



also references for the one-channel case |2^). But there exist SUSYQ transforma- 
tions which are non-equivalent to the IP ones. For example, one can accomplish the 
double-SUSYQ procedure at different energies so that the first step is made at a 
bound state energy while the second one is performed at another energy, see, e.g.. 



for the one-channel case the work p4| . Then a chosen bound state energy is shifted 
and the spectral weights for all the energy levels are changed. As infinite number of 
spectral parameters are involved, to fulfill such a transformation is infeasible in the 
framework of the IP. Thus, in some cases, the SUSYQ formalism and one-parametric 
IP lead to the same results, whereas in some cases, SUSYQ transformations are not 
equivalent to finite-parametric models of the IP. 

The interaction matrices V"q,^(x) can also correspond to non-local potentials 
V{x, ^, ^') of the related multi-dimensional problem. So, wide classes of exactly 
solvable models mentioned above open new possibilities to investigate systems with 
non-local forces which are still very poorly understood. 

3 Uncoupled channels 

Let us start the consideration of the multichannel peculiarities with the limiting 
case of disconnected equations (Q) when Va^p = 0. The partial channel transformed 
potentials and the corresponding wave functions can be unified in the matrix, e.g., 
2x2 form 



n^) = 'T' . (13) 
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(14) 



For uncoupled channels it is especially clearly seen that the spectrum consists of 
several branches (each for its channel). By altering a chosen threshold value ea, 
one can shift the ath partial branch with respect to the others. Some bound state 
energy level of one spectral branch can occur in the continuous spectrum of another 
channel. It is the simplest case of a bound state embedded into the continuum 
(BSEC). In the general case of coupled channels, the spectrum becomes unified and 
its branches are not so clearly separated. 

It is interesting that one can join the initially independent branches at energy 
En through creating at this point a common bound state (or BSEC if En is above 
at least one of the thresholds). This is attained in the IP formalism when both 
components Can of the SWV for E = En are chosen to be non-zero. It results in 
some channel coupling Vi2ix) according to, e.g., Eqs. (36)-(38), section 4.7. At 
other energies, spectral characteristics remain unchanged and correspond to initial 
separated channels despite arbitrarily strong coupling. And only at the chosen 
energy En the spectra are mixed completely. For more details on the SWV control, 
see section 4.3. 

The system of uncoupled channels is a convenient bridge linking one- and multi- 
channel theories. Here the one-channel 'bricks and building blocks' p, ^ [Tl|, |12|, ^ 
appear in their pure form. 

As the first example, we shall consider the two-channel model with infinite rect- 
angular potential wells Vii(x), V22{x) having purely discrete partial spectra. The 
decrease in the spectral weight parameter Cu for a bound state in the first channel 
is demonstrated in Fig.l. This makes the chosen normalized partial wave smaller at 
the origin and bigger on the right due to conservation of the norm. It means shifting 
the chosen partial state in the configuration space (along x) to the right. The cor- 
responding potential transformation A Vii (x) has a barrier on the left which pushes 
the wave to the right. But this barrier alone would shift upwards all the energy lev- 
els of the first channel and violate the restriction that all other spectral parameters 
except for Cu remain unchanged. To compensate this undesirable influence of the 
barrier, the additional well is introduced on the right. 

A general rule-I for the spatial shift of an arbitrary partial state in any poten- 
tial was found §. By decreasing (increasing) Can the nth partial state 
is shifted to the right (to the left) by n analogous Al^Q(x)-blocks: 'bar- 
rier-|-well' ('well-|-barrier') one for each bump of the wave function. Re- 
turning to the initial system obviously requires the inversion of the perturbation 
block: the first transformation (barrier+well) combined with the second (the reverse 
'well-|-barrier') cancel one another and restore the initial rectangular potentials. 

The increase (decrease) in the spectral parameter Cu at one boundary of the 
allowed interval of motion is equivalent to decrease (increase) in the adjoint spectral 
parameter Mu defined at another boundary (or for x oo). The same is not valid 
when an interchannel coupling is switched on (see below). 

Shifting a chosen energy level over the ath branch of the spectrum is another 
group of the complete set of elementary transformations. The case of potential 
symmetry with respect to the center of the interaction interval is especially simple 
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[|, H, |2J]. Here, keeping the spectral weights unaltered is superfluous because noth- 
ing but energy levels form a complete set of spectral parameters for an uncoupled 
channel. 

The qualitative rule-II. Upward (downward) shifts of the energy level En 
require barriers (wells) in the potential AV^o(x) at the regions of the most 
sensitivity of the nth chosen state: near the extremum of each anti-knot. 
Compensating wells (barriers) at positions of knots are necessary to keep 
other energy levels at the same positions. It is illustrated in Fig.l for the 
lowest bound state in the second channel. 

Let us give a notion of analytic expressions for these transformations. The for- 
mulae for the transformed potential and the wave function of the lowest bound state 
related to the case of changing Cu are : 

(.) - 2±. . (15) 

''^ l + (CiV C„-l)/„"=*i(!/M!/ 

and 

i + (C?i/c„ -i)SiAy)iy 

o o 

It is remarkable that the new system is constructed of the initial (x) and Vii {x) 
as building material. Analogous simple formulae for the one-channel energy shift 
are given in |2^ . 



The same rule for spatial shifts is valid in the case when there exists continuous 
spectrum. If the external potential wall is of finite height, the last well of the nth 
block can be pressed out of the initial potential (carrier potential wells of the soliton- 
like form - the fundamental quantum brick). It is shown for the second channel in 
Fig.2. 

In the limit Can — ^ or M^n oo, the bound state is carried away to infinity 
inside the refiectionless well. It is equivalent to the effective removal of the chosen 
energy level from the initial spectrum . 

A scattering state wave can be considered as a bound state in the infinite broad 
potential well. If we imagine it to be distributed over the infinite interval with norm 

o 

1, it has an infinite small 'spectral weight' CaE- By making CaE finite, this state 
can be raked up to the origin and becomes BSEC. But there is an infinite number 
of bumps in the initial state. So we need an infinite set of the potential well-barrier 
blocks, as shown in Fig.2. It is the limiting case of the above-mentioned rule-I. The 
BSEC in one channel can coexist with a scattering partial component in another 
channel at the same energy. It is a non-trivial fact that this can be realized for 
strong coupling of channels as well, see below. 

It is remarkable that energy levels can also be shifted in the "imaginary direction" 
[ |r^ : E E + iV. The above stated qualitative rule is valid for Im^Vaa{,x). So 
quadratically integrable decaying states (with the exponential factor exp(— Ft)) are 
new convenient exact models. By analogy with the optical model it is an effective 
account of coupling with virtual channels not exhibited explicitly. 
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4 Coupled channels 



As in the one-channel case, we can perform here a complete set of utmost elemen- 
tary transformations. Any chosen bound state energy level can be shifted and any 
partial channel component Can of the SWV can be changed using the closed analytic 



expressions, see ||13| , P5|1 and references therein. We shall consider some qualitative 
properties of multichannel solutions after some preliminary remarks. 

Switching on non-diagonal elements of the interaction matrix in leads to a 
wave exchange between the channels which resembles a fluid motion in communicat- 
ing vessels. The conservation of a partial channel flux is usually violated and only 
the total current is conserved. 

Below are some simplest examples illustrating the specific influence of channel 
coupling. If we introduce a constant coupling V12 = V21 = W = const into two free 
equivalent equations of motion with initial solutions ~ sin(/cx), k = \/E 

-^'2'(a;) + W^i{x) = ^^2(3;), (17) 

we get new solutions with frequency that depends on the boundary conditions. The 
partial waves coincide = ^1/2(2;) k = \/E — W for identical partial boundary 

conditions. But if the partial boundary conditions are chosen with different signs 
the solutions = — \E'2(x) ~ sm{kx) have another wave number k = \/E + W. 

Let us consider a system on the whole line with equal thresholds and rectangular 
interaction matrix elements Vii(x) = V22{x) = —Vi2{x) = V, < x < a (Va/six) = 
0, X < 0; a < x) 

- ^"(x) + V^i{x) - V^2{x) = E^i{x) 

-^2{x) + V'^2{x) - V^i{x) = E^2{x). (18) 

In this case, we have the absolute transparency for identical boundary conditions 
in both the channels \l/i(a;) = '^2{x)- Moreover, there is free (!) motion due to 
complete cancellation of interactions. But if we choose different signs of the solutions 
\l/i(x) = — \l/2(x), there is no such cancellation, which results in strong reflection of 
waves. 

The next example demonstrates the influence of channel coupling on energy 
level spacing. Switching on the constant coupling W between two channels with 
equivalent thresholds and infinite rectangular potential wells Vu = V22 results in 

o 

splitting initially degenerated levels En =En 

Now let us consider the multichannel spectral control by using the IP approach. 



4.1 Variation of SWV components 

o 

Let us start with an arbitrary initial potential matrix Va/3 (x) for the system with 

o o 

bound states. We shall vary SVW Can^ Can and energy value En^ En of the 
chosen nth bound state. The relevant expressions for the transformed potential 
matrix and solutions are derived, e.g., within the framework of the multichannel GL 
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formalism. This is the case that corresponds to the exactly solvable models of IP 
discussed in section 2. Similar formulae are in the Marchenko approach as well. 
The resulting potential matrix has the form 

Vaf,(x) ^VaP {X) - 2^{f(x)P(x)-lt(x)U, (19) 

where 



/ o o o o 

^(^\- ( E/3 Cpn $1/3 {X, En) E/3 C/3n$l^ {x, En) 
^ \X ) — I o o 

E/3 Cfjn $2/3 {X, En) E/3 C/3n$2^ {x, En) 



(20) 



and 



Y(2;) ^ I E/3C/3n $1/3 (a:, -En) E/3 C'/3ri $2/3 (a:, -E„) i ^^^-^ 
- E/3 C/3n$l/3 {X, En) - E/3 C/3n$2/3 {x, En) 



Pu{x) = 1 + riYCpn $1/3 {y, En)YCpn $1/3 (l/, K) 
"'^ /3 /3 

+ C/3n $2/3 (y, En) Pn $2/3 (l/, En)]rfy 

/3 P 

Pl2{x) = -P2l(2;) = TE $1/3 (?/, K) E &n$l/3 (?/, In) 

/3 /3 
_ — ^ o ^ — ^ o o o 

+ E ^P-^ ^n) E C'/3n$2^ (?/, -En)]c?y 

P P 

f'^ ^ o o o ^ o o o 

P-22{x) = 1 - / [Z^ C/3n$l/3 (y, -En) ^ C'/3n$l/3 (?/, En) 



P 

o o 



+ E C/3n$2/3 (y, En) E C/3n$2/3 (|/, En)]c?y. (22) 

The transformed regular solution at arbitrary energy E is 

^^p{x, E) ^hap {x, E) - r Y.{\{x)P{x)-^T{y)}^^ h^p {y, E)dy. (23) 

The solution \E'Q,(a;, En) = E/3 C'/3n$a/3 (a;, En) = {T{x)P{x)^^}ai represents a nor- 
malized wave function of the pushed nth bound state 



/ Yl\.'^^{x,En)]'dx^l. 

Jo ^ 



Computer visualization of these formulae shows curious peculiarities of the corre- 
sponding transformations of the potential matrix and wave functions some of which 
are presented here. 

A gradual increeise in only one peirtial channel component of the spectral 
weight vector Can or Man for the nth bound state results in progressive 
concentration of waves in the ctth channel by their transition from all 
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other partial channel components j3 ^ a oi the whole nth wave function. So 
the waves are gathered from botli the configurational and channel spaces while the 
total norm is conserved. In the limiting case Can — > oo or Mq,„ — > cxd, the whole 
wave function is pressed into the origin in the a channel or is shifted to infinity. 
It means that the chosen energy level effectively disappears from the spectrum. A 
typical behavior of Vaf3{x) and the wave function is shown in Fig. 3 when Mn oo. 
The wave function is gradually moved away to infinity by the refiectionless soliton- 
like 'carrier' potential well. This process is accompanied by the new effect: the 
concentration of all waves in the chosen (first) channel at the expense of the others 
(second one). This is true for both equal and different channel thresholds. What 
is also amazing is that in spite of strong coupling of channels, all other channels 
are almost completely emptied. The 'carrier' interaction well is surprisingly similar 
to the one-channel case: we meet again with the soliton-like potential as a simple 
"brick" of quantum reconstruction. The discovery of this remarkable effect deepens 
our understanding of complex multichannel systems. So in comparison with the 
one-channel case, there is a strong magnification of the response of the system to 
increasing partial components Can or Man- 

On the contrary, the gradual decreasing to zero of only one partial chan- 
nel component of spectral weight vector Can or Man does not completely 
empty the ath channel, see Fig. 15 in |jl2[ as a typical example of such a situa- 
tion. The ath channel wave is only partly shifted. Some part of it is pushed out into 
other channels. In the one-channel case, the nth state disappears if C„ — >■ 0. In the 
multi-channel case, the zero value of the partial function and its derivative at one 
point does not mean that this wave must disappear. It returns to this channel at 
other points due to the wave exchange between channels. Only if all the components 
Can become zero, the bound state is completely removed. Thus, in comparison with 
the one-channel case, we have the weakening of the response of the system to the 
decrease in the partial spectral weight parameter. Furthermore, the equivalence of 
increasing (decreasing) in Can and decreasing (increasing) in Man being valid in the 
one-channel case is violated when the channels are coupled. 

The IP and SUSYQ formalism allow one not only to remove energy levels but 
also to create them at given positions. With the initial free motion system this 
creation leads to refiectionless interaction matrices. 



4.2 Transparency 

It was once found that the Coulomb barrier between nuclei with large Z is much 
more penetrable than was expected. Before our investigations on IP, we found the 
phenomena of intensified barrier transparency and even supertransparency in the 
multi-channel approach to quantum many-body systems, see e.g. 0]. 

Additional understanding of the phenomena can be achieved via consideration 
of the limiting model of total transparency. The isospectral transformation of a free 



motion system into the refiectionless one with a bound state gives [26] 

d MaMfsexp[—{Ka + Kjs) x] 

-7 2^^^P(~2/t7 



Va^ix) = 2- ° ' ^ " \ Ka = ' E,, (24) 
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and 

^a{x,Eb) = jf^ -. (25) 



Here corresponds to the created bound state. In |2^, the transparent interac- 
tion matrices in the case of different thresholds were derived and shown. For equal 
thresholds, the matrix elements of Vaa{x) have a simple soliton-like form. For differ- 
ent thresholds, there occurs an repulsion in Vaa{x) which is unexpectedly necessary 
for the complete transparency. The waves reflected disappear via destructive inter- 
ference with the backward waves decaying from other channels. It is a non-trivial 
multichannel analog of the one-channel soliton-like potential well. 

In (^), the first channel component of the bound state wave functions decreases 
"unnaturally" differently in the directions x ±00. It may seem that the channels 
become "disconnected" at large values, and partial waves must have a natural 
asymptotic decrease expi^^y/E — eQ,|a;|), x — >• ±00 for Va^{x) rapidly decreasing in 
both the directions. But even a weak coupling Vaf3{x) — >• as x — > cx) can suck 
out the remainder waves from some channels into other ones violating the standard 
asymptotic behavior. Instead of the "naturally expected" behavior 

'ifa{x,Eb) — >■ exp{—Kax), X ^ oo (26) 



we have 



and 



^i(a;,Eb) ^ exp[(-fi;i + 2^2)0;], x ^ -oo (27) 



^2(2:, -Efc) exp{K2x), X —00. (28) 

So, there is a phenomenon of partial inversion of the " degree of closeness" of different 
channels. 

This can also be illustrated by reducing the two-channel system (0) to one equa- 
tion for the partial channel wave function \E'i(x, i?f,). This is done by substituting 



the explicit expression for \I'2(x,-E;,) in terms of '^i{x,Eh) (see Eq. (|25|) ) 



Mo 

^2{x, Eb) = ^ exp{(-K2 + Ki)x}^i{x, Eb). (29) 

into the first equation in (0) with the potential matrix (p^. As a result, we get the 
one-channel Schrodinger equation for \E'i(x, i?b) with the effective potential 

V[x) = Mfexp[-2(/ti) x] ^ ^ 

l + E^I?-exp{-2K^x) 



,d M.M.exp|-(.. + ^^^ ^^^^^^ 

l + E7^exp(-2/€^a;) ^1 



As X ^ —00, it becomes a non-zero constant 4k2(k2 — /ti) > whereas ^22(2;) 
decreases as exp[(2/t2 — 2fi;i)x], which leads to the unnatural inversion of closeness 
of channels (on the left). This illustrates the fact that the rate of the asymptotic 
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decrease in the partial channel wave function can be strengthened even by exponen- 
tially decreasing coupling V^^. The interaction ( pO]) is energy dependent. 

In [|l6l, we have also found transparent interaction matrices without bound states 
by employing the multichannel SUSYQ formalism which gives a broader class of 
exactly solvable models. These quantum systems not previously known have no one- 
channel analogs. They play a significant role in creating a bound state energy level 
close to the existing one in the initial spectrum (see the next section and the left 
part of Fig. 5). 

One of the multichannel peculiarities found by us is the right-left symmetry 
violation [§, ^ for the reflection of a complex particle by a non-symmetrical po- 
tential barrier. It is more penetrable in one direction for an incident particle in the 
ground state. And in the opposite direction (in two-channel approximation) it is 
more transparent in the exited state. 

Let us now consider the model demonstrating shifts of energy levels up to their 
degeneracy. 



4.3 Degeneracy of energy levels 

In the one channel case, two bound states are forbidden to have the same energy. 
We have discovered the phenomenon of "effective annihilation" of degenerating 



states. These states must both be equal, as solutions of the same equation with the 
same boundary conditions, and simultaneously orthogonal as different states. This 
contradiction is resolved in the following way. There is a separation of parts of the 
waves by the non-penetrable barrier (infinitely high or broad) under which \E'q, — > 0. 
So the wave functions become identical up to the sign in the limit of zero distance 
between the levels En — En+i — > 0. But in the case of M channels, M degenerate 
states can coexist provided that their SWV are hnearly independent. 

Let us use the IP formulae for creation of two energy levels Eh and E'^ which 
are very close to one another (including the case of degeneracy). As an initial two 

o 

channel system on the whole line we choose the free motion with Vafs (x) = 0. These 
formulae are similar to those in section 4.1: 

Vo^^ix) = 2^{f(x)P(x)-it(x)W, (31) 



where 



^ ' \M2exp{-K2x) M^expi-K'^x) )' ^ ' 



Tix) 



Ml exp(— Kix) 
M[ exp(— k'^x) 



and 



M2 exp 
M^expi 




(33) 



Piiix) = 1 + 



E 



M2 



exp( 
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P,,{x) = 1 + Y: ^ exp(-2<a:). (34) 



Here = V^a — Eb, i^'a = y^a — E'^ and M^, are the SWV of the states under 
creation with the energies Ef, and E'^^, respectively. Let us combine normahzed 
vector-column wave functions at the energies Eb and E^ to form the matrix \l/(x) 
for compactification of the notation 

^{x) = f{x)P{x)-\ (35) 

When two levels Ei, and E'j^ with linearly dependent SWV come closer and closer 
to each other, the features of the effective annihilation become apparent. The wave 
functions are split and removed (either their parts or the whole states), as in the 
one-channel annihilation case. The same happens if SWV of the degenerated states 
become more hnearly dependent (see Fig. 4). Rudiments of spatial separation (prepa- 
ration to the effective annihilation) can be observed even without changing the en- 
ergy distance between the levels. As in the one-channel case {^^, Fig. 13), the 
states with linearly dependent SWV even corresponding to different energy levels 
will resist the concentration in the narrow spatial region. There was found an ex- 
change of knots between partial channel wave functions with flips of some of their 
bumps (§, Fig. 14). The annihilation phenomenon and the recoil of all states with 
quantum number n ^ m from the origin when the chosen mth state is concentrated 
to X = p|, |25[ have much in common. They are connected with the resistance of 



system of different states to their concentration in phase space. 

There also occurs 'splitting' of BSEC when Eibsec E2bsec (as has been 
shown in [ig. Fig. 22). 



Unexpectedly, when two levels even with linearly independent SWV come close 
to each other, this is often followed by special evolution of the interaction matrix 
Vai3{x). In the case of unequal thresholds, there occurs a spatial separation of some 
block of Vapix), see Fig. 5. This block is transparent and goes to infinity in the limit 
of the level degeneracy resembling one-channel annihilation This separated 

block, in contrast with the annihilation case, does not carry away any part of the 
bound states. All bound states remain completely inside the main part of Vajs^x). 
We have already mentioned this peculiar interaction matrix in the previous section. 
An explanation of this new effect is still an open problem. 

It is interesting that, in the multichannel case, there is the possibility of moving 
an energy level through the other level positions without the appearance of singu- 
larities in the interaction matrix. Such an operation is feasible if the SWV's of the 
shifted state and the crossed states are linearly independent. The crossed states 
must be less then M-fold degenerate where M denotes the total number of channels. 

4.4 Bound states embedded into the continuum 

BSEC being a wonderful phenomenon in the one-channel case, has a diversity of 
new multichannel aspects. A possibility appears to control the asymptotic decrease 
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of BSEC: making BSEC short- or long-range. The examples of exponential BSEC's 
falloff Va/3{x) and -Ebsec) for systems with BSEC below the threshold of the 

most closed channel were shown in [^]. Later, we have discovered that in special cases 
BSEC can decrease like ~ 1/x when its energy is between the threshold ones. It 



can be explained as follows. The regular solutions increase exponentially in closed 
channels as x — oo. The physical scattering solutions i?) are constructed 

as a linear combination of the columns of the matrix of regular solutions with the 
coefficients providing exponential decrease in the closed channels. Let us choose the 
SWV components for BSEC proportional to these coefficients for solutions in the 
initial system (with non-vanishing channel coupling) where BSEC is to be created 
at Ebsec- Then the resultant Va/six) and "^aix, Ebsec) have ~ 1/x behavior. It is 
caused by the linear increase of the integral 



^ rx o 2 



in the denominator of the formulae for Vaj3{x) and "^aix, Ebsec) which are direct 
GL analogs of Eqs. (36) - (38), section 4.7. In this expression exponentially growing 
terms are cancelled by the aforementioned choice of SWV components Ca- Any 
violation of the ratio between Ca leads to the exponential growth in the denominator. 
This results in exponential falloff of Va/3{x) and "if a{x, Ebsec)- 



As in the one-channel case |^, a BSEC solutions above the thresholds, despite 
the different frequencies of oscillations of partial wave components, have strict con- 
formity between knots of functions and diagonal interaction matrix elements. For 
each partial function bump there is a relevant potential well-barrier block, see Fig.l 
and |2|. 



We have considered in section 4.1 pumping waves into a chosen channel by in- 
creasing the corresponding component of SWV. Analogous effect can also be ob- 
served for BSEC. The BSEC can be obtained even by taking a single partial ath 
component of SWV to be non-zero. Then the interchannel coupling provides the 
transformation of other partial scattering waves into BSEC. 

4.5 Phenomena of coexistence of states with 'incompatible' 
properties 

The existence of several physically allowed linearly independent solutions of Eq. (0) 
make it possible, unlike the one-channel case, to combine at the same energy different 



peculiarities of wave dynamics. For example, in |[T6| it was revealed that a bound 
state can coexist with a scattering state even without a strict space separation of 
propagating and bound waves. The qualitative arguments to resolve this "paradox" 
are the following. In the M-channel case, there exist M linearly independent regular 
solutions of (0) for the same energy point. Above all thresholds these solutions are 
scattering states. In the IP and SUSY approach one of the scattering states can 
be transformed into a BSEC aX E = E^ > ef, i = 1,2,...M ). Then there will 
remain M — 1 scattering state solutions. As a result, it is impossible to construct 
scattering solutions at BSEC energy with arbitrary asymptotic conditions, e.g., with 
the incident wave in one channel only. 
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By analogy with the effect considered above, there can exist scattering solutions 
with different resonance decay widths at the same energy. We can construct such 
an exactly solvable model. Let us start from two uncoupled channels with different 
thresholds and resonances at the same energy with unequal widths ri,r2. We can 
create a bound state with non-zero SWV components in each channel by using the 
GL analogs of Eqs.(36)-(38). There arises a strong coupling V^ix) between channels. 
For incident waves in the ath channel only, we shall get resonance scattering with 
the width Fq, despite the intensive exchange of waves between the channels inside 
the interaction region. Outgoing waves survive in the entrance channel only. After 
the mixing of channel waves in region where Vi2{x) ^ 0, they return to the entrance 
channel. It is a common opinion that compound states (resonances) in complex 
systems "forget" how they were generated. Here an extremely simple example has 
been demonstrated of direct dependence of the life-time of a quasi-bound state on 
boundary conditions. So in the general case the analogous phenomenon is also 
possible. 

In section 4.7, the coexistence of transparency and strong reflection is shown. 



4.6 M+1 spectra 

The systems with the pure discrete spectrum in the one-channel case are uniquely 
determined by the set of C^} or by two spectra {EniEm)- Here the index 
m is related to eigenvalues with different boundary conditions at one end point 
of an interval (e.g., at the origin \l/^(0) = instead of \l/a(0) = 0) 0, p. 401. 
Daskalov has suggested exactly solvable models corresponding to variations of 
any eigenvalues from two spectra -Em}- This enlarges the possibilities to control 
quantum systems. In the M-channel case, the set of eigenvalues and the sets of M 
SWV components {E^, Can} form the complete (M+l)-fold set of parameters. They 
are equivalent to M-l-1 spectra with different boundary conditions [§]. Exact models 
for M+1 spectra variations are now under construction. 



4.7 Resonance tunneling 

Almost non-penetrable barriers could be combined in a system transparent at dis- 
crete energy values, see e.g. This phenomenon in the one-channel case has 
already found numerous applications, see and references therein. But we have 



not seen any papers on multichannel resonance tunneling. The waves accumulated 
between the barriers in the Vapix) decay in both directions. The decaying waves 
going backward can have the opposite phases and the same amplitudes as the re- 
flected waves in all open channels. In this case reflected and decaying waves cancel 
one another totally, which results in resonance tunneling. It appears that the multi- 
channel analog of this phenomenon is restricted by additional conditions. We have 
found that it can be achieved only for special proportions of amplitudes of incident 
waves in different channels at discrete resonance energies Eres- 

So, there can coexist at the same energy solutions without reflection at all and 
with strong reflection depending on boundary conditions. The relevant example (ex- 
actly solvable IP model) can be constructed of the two initially uncoupled channels 

o o 

{Vi2 (x) = 0). We can choose the potential Vu (x) in the first channel having 
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"one-channel" resonance tunneling at E = Ej-^s and in the second channel the po- 

o 

tential V22 (x) being weakly penetrable at the same energy. Let us now create 
a bound state at E = E^ < common for two channels by the IP transforma- 

o 

tion Vai3 (x) — * Vai3{x) that mixes intensively partial channel waves (i.e., leads to 
^12(3;) 7^ 0). Then the non-trivial two-channel system occurs for which the waves 
incident in the first channel (i.e., ipaix) 5ia exp(— z/cix), x 00) propagate with- 
out reflection at E = Ej-es whereas the waves incident in the second channel (i.e., 
'ipa{x) exp{—ik2x), X oo) are strongly reflected. This is possible because the 

IP transformations keep the continuous spectrum characteristics unchanged. The 
corresponding analytic expressions have the form : 

o o 
/ X ° f \ , Fafi' {x,Eb)Ma,Mf3, F (3 a' jx, Eb) 
Vafsix) =Vaf3 [X) + 2- ^ , (36) 

1 + C E-y Ea'f^' F^p' (y, Eb)Ma'Mp, Fya' {v, Eb)dy 



where = \/ea — E^ 



Fap{x, E) =FaP (X, E) - 



E -^-/^'(;-^^)^"'^/^' F,'a'{x'^Eb) ^^,^(,.^)^,. (37) 

y 1 + /r E-, Y.a'p' Fjfs' (y, Eb)M^,Mp, F^a' {y, Eb)dy 

where Fapix^E) is the Jost solution at arbitrary energy E obeying the asymptotic 
condition Fap^x^E) Sap exp(— zv^E' — ea x) as x ^ oo. The created bound state 
has the form 

E,) ^ ^f,M,F„,(..E.) ^ 

1 + XT E-y Ea'/3' F^p' (y, Eb)M^'Mp, Fja' {y, Eh)dy 



We have found in ^3] that two coupled particles may exhibit resonance tun- 
neling properties while undergoing transition through only one potential barrier. 
This problem can be reduced to multichannel equations with two barriers in the 
interaction matrix. So, our consideration above applies in this case as well. 



4.8 Periodic interaction matrices 

As a particular limiting case of resonance multichannel tunnelling one can consider 
the wave motion in periodic fields. In this case the discrete resonance tunneling 
points for two-barrier system broaden into allowed continuous spectrum bands. 

To develop the theory of excitations of crystals, there can be useful exact solu- 
tions for quasi-particles with the inner (multichannel) degrees of freedom. 

Let us consider the two-channel system. For two uncoupled channels with pe- 
riodic Vaa{x) there are two independent spectral branches with energy bands of 
allowed and forbidden motion with quasi-momenta Ki 2. For example, in the case 
of Dirac's comb potentials Vaa{x) = I]^_oo Va5{x + na), each of the channels we 
have 

sm{kaa)Va f, \ u rQn^ 

cos(AQa) = h cos[kaa), ka = y E — e^. (39) 
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One could suspect that switching on a weak couphng would hinder wave propagation 
at the energies belonging to a forbidden zone in one channel and to an allowed zone in 
another one. Indeed, it may seem that the exponential increase in the wave function 
component in the forbidden zone would prevail and cause physically unacceptable 
asymptotical divergence of solutions. But really it appears that switching on the 
coupling changes slightly the allowed zones of the uncoupled channels. So just the 
regime of oscillations, characteristic to the allowed zones, dominates. It results in 
damping the exponential increase corresponding to forbidden zones. In fact, there 
are two types of two-channel generalized Bloch's waves 

^^Bai^) = exp(±z/C(^'2)a)^gf - a) (40) 

with the quasi-momenta /C^^'^^(£') common for both the channels. We can find 
such quasi-momenta by using, e.g., the exact model with the following periodic 
^-interaction matrix : 

oo 

Va/3ix)= Va(35{x + na). (41) 

n=— oo 

With such a choice of K,/3(x) we have 

cos(/C*-^'^''a) = -{cos(A'ia) + cos(i^'2o) 

L 77^ ^ 77^ sin(A;ia) sin(A;2a) _ , , 

±J[cosiK,a) - cosiK^a)]^ + ^ ^ ^ }, (42) 

where W = V\2 = V21, Va = Vaa, and the expressions for cos(-ft'i_2a) are given by 
(PUD- So the condition \ cos{}C^^''^^ a)\ < 1 serves for the criterion of searching the 
allowed zones. It is shown in Fig.6 that the quasi-momenta /C^^'^^ are close to ones 
for the uncoupled channels (for W = 1). 

Another peculiarity is that we can create a gap in the given allowed spectral 
zone constructing a special multichannel periodic potential. In fact, let us consider 
at first wave motion on a finite interval (period) [0, a] with homogeneous boundary 

o 

conditions. Further, we transform a constant interaction matrix Vai3 defined on this 

o 

interval to Va/3{x) so that a chosen bound state at the energy En related to Va/3 (x) 
is raked up to the right boundary by scalar decrease in the SWV Can- Next, let us 
continue the resulting matrix block to the whole axis with the step a 

V^'/{x + la) =Vap{x), / = 0,±1,±2..., 0<a;<a. 

For this periodic potential, Bloch's solution at the energy En can be constructed 
in each Ith interval [la, {I + l)a] from the solutions obtained via the /a-transfer of 
"i/aix, En) over the coordinate line: \l/^^-*(a; + la, En) = \l/a(x,-E„), < x < a. We 
have the following equalities: ^„(0, En) = "ifaia. En) = and ^i(0, E„) 7*2(0, En) = 
\l/'i(a, En)/'^2{ci, En). This guarantees the possibility of smooth continuation of the 
wave function under the construction *^^^(x, E'n) (—00 < x < 00) at the points 
la provided that the module of the solution ^^'^^^x. En) for the (/ + l)th sub- 
sequent sector is each time multiplied by a scalar (i.e., independent of a) factor 
e = %ia. En) /%iO, En) > 1 (« = 1,2) 

*P^'^(x + /a,E„) = (-l)'e'*„(x,E„), / = 0,±1,±2..., 0<x<a. 
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Eventually, this leads to an exponential swinging divergence of En) as x ^ oo 

at the energy En- It means that this energy belongs to a forbidden zone created 
with the above procedure. 

We can transform the periodical multichannel system by creating bound states 
inside the allowed or forbidden zones as, e.g., was shown in It is also interesting 
that using the SUSYQ transformations of (41) we have an additional flip effect of 



all the delta-peaks, as in the one-channel case [34]. In fact, after the single SUSYQ 



transformation (see, e.g. [|T6]) the potential matrix has the form 

V{x) = Vo{x) - 2^[^'o(x)^o(x)-^] = -K)(x) + 28 + 2[^[,(x)^o(x)-T, (43) 

where Vo{x) is given by (0); "^oix) represents the corresponding matrix- valued 
solution at the factorization energy S. The term —Vq{x) with the minus sign occurs 
in the last equality, which means the flip of all the original (5-peaks which cannot be 
compensated for all x by other finite terms. 



5 Conclusion 

The quintessence of the quantum mechanics is the special relationship of observables 
S and interactions V for wave systems. The universal simple rules for V elementary 
transformations caused by elementary S variations in the one-channel case were 
established earlier p, . 



For a long time the multichannel theory has been a 'black box' connecting input 
and output data through cumbersome computer calculations. The examples pre- 
sented in this paper help one to look into intimate physical details of different new 
effects inherent in coupled Schrodinger equations. 

The considered elements of the multichannel quantum design give some notion 
about elementary transformations from their complete set for systems with the com- 
plicated inner structure. No doubt, the multichannel theory must be richer and more 
wonderful than in the one-channel case. 

From a historical point of view, it can be worth mentioning that the IP for multi- 
channel systems with the same thresholds was considered in the book by Agranovich 
and Marchenko on multichannel IP (see references in [^). See also the interesting 
papers [|15|, |l^ |18], |3^ and papers by Cox (see references in [E|). 
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FIGURE CAPTIONS 



Fig.l. Two uncoupled channels with the initial infinite rectangular potential 
wells Vaa{x) transformed in the following way. The lowest bound state in the first 
channel is shifted over the coordinate axis x. The lowest energy level in the second 
channel is shifted over the E axis. The ground state of the first channel spectral 
branch was shifted to the right. All the energy levels in this branch and all other 
spectral weights C„^i remain unperturbed. To decrease only one spectral weight 
ip'i{x = 0) of the ground state of the first channel, a potential perturbation block 
AVi(x) (barrier + well) is needed. The lowest energy level of the second channel 
branch was lifted upward without shifting other levels of the second branch. 

Fig. 2. A scattering state is transformed into a bound state at the origin in the 
first uncoupled channel by increasing its spectral weight C{Ebsec)- 

Soliton-like well in the second channel serves as a carrier for the chosen bound 
state moving out of the original well when the partial channel spectral weight vector 
component decreases. 

Fig. 3. Increasing first channel parameter (Mn = 10 Mii,M2i =M2i) for the 

o o 

ground state. Both the initial wave functions {x,Ei) and \[/2 {x-,Ei) are concen- 
trated after the transformation in the first channel wave component '^i{x,Ei) and 
pooled to the right inside the separate soliton-like potential well Vu {x) . The initial 

o o 

potential matrix is defined on the half-line and has the form Vu {x) =V22 (x) = 

o o 

—5, Vi2 (x) =V2i (x) = 0.3, < X < TT, Va/3{x) = 0, X > 71, the thresholds are 
ei = 0, €2 = 1. The remaining part of \l/2(a;,i?i) in the second channel is so small 
that is not visible. 

Fig. 4. Two degenerated bound states created at Eb = —0.5 become 'more' 
linearly dependent. This results in splitting of a) the interaction matrix and b) 
wave functions into two parts. The initial system corresponds to free motion along 
the whole line, a) Two blocks of Vap{x) for almost equal SWV: Mi = M2 = 1; M[ = 
1, M2 = 1.01. The left part of Va/six), a soliton-like well 'moves' to —00 in the limit 
M2 — M2 as is schematically shown by dashed lines, b) This well carries also the 
left part of wave functions to —00. 

Fig. 5. a) A potential block (left part of the interaction matrix) without a bound 
state is separated and shifted to the left when two closely spaced bound states are 
created : = 0.5; Mi = 0, M2 = 1 Eb' = 0.501; M[ = 1, M^ = 0.1. This block 
'moves' to — cx) in the limit Eb — ^ Eb' as is schematically shown by dashed lines, b) 
Wave functions remain almost unchanged. They don't "live" inside this potential 
block, which is almost transparent if considered separately of other part of potential 
matrix, see and Fig. 8 there. 

Fig. 6. a) The functions cos(/C*^^'^^7r) (solid fine) specifying the location of spectral 
zones for the two-channel system are shown. Allowed and forbidden zones are de- 
termined by I cos(/C'^^'^-'7r)| < 1 and | cos(/C''^'^''7r)| > 1, respectively. The parameters 
are Vi = 6; V2 = 5; = 1; a = vr; ei = 0; 62 = 1. In the case of uncoupled channels, 
the corresponding curves cos(i^i^27r) for separate channels are shifted from one an- 
other due to the threshold difference and have crossings shown by the dashed lines. 
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Switching on the coupling transforms the crossings (dashed hnes) into the quasi- 
crossings so that there are the upper and the lower solid lines for two branches of 
band spectra, b) Solid line intervals are forbidden spectral zones of the coupled equa- 
tions (their allowed zones are outside these intervals) . The sohd and the dashed hne 
intervals together represent the forbidden zones for any uncoupled equation. Only 
blank intervals between vertical dashes correspond to the common allowed zones for 
both uncoupled equations. 



22 




Fig.l 




Fig.2 




Fig.3 




Fig.4 




Fig.6 



